The inversion-asymmetry of pion emission source along the 
outward direction in relativistic heavy ion collisions 



o ■ 
o ■ 

(N : 

s . 

<: 

m : 
> 

o : 
\o . 
o ■ 

00 1 

o : 
\o ■ 
o ■ 



o 

=5 



X 



Liu Lianshou* Shi Shusu and Du Jiaxin 1 

Institute of Particle Physics, Huazhong Normal University, Wuhan 430079, China 

The inversion-asymmetry of the pion emission source in relativistic heavy ion collision under 
the Bertsch-Pratt convention is discussed and explicitly exhibited by a Monte Carlo model. The 
Gaussian source function popularly used in the HBT analysis of relativistic heavy ion collisions is 
invalid in this case. An inversion-asymmetric source function is suggested. A method for extracting 
the inversion-asymmetry degree of the source together with the source size from experimental data 
is proposed. 
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The analysis of two-particle, in particular two-pion, 
correlation provides a powerful tool for determining the 
spatio-temporal characteristics of relativistic heavy ion 
collisions yj. It is usually referred to as pion interfer- 
ometry or HBT analysis, following the pioneer work of 
Hanbury-Brown and Twiss |2(. Unfortunately, there is 
no one-to-one correspondence between the two-particle 
correlation function in momentum space and the source 
distribution in configuration space. How to extract as 
much information as possible about the source distribu- 
tion from the experimentally measured two-particle cor- 
relation data is the main challenge for the HBT analysis. 

Two-particle correlation function can be expressed as 



C(P,q) = l+ / d 3 rS P (r)|0,(r)| s 



(1) 



where P = p\ + P2 is the total momentum of the pair 
and q = pi — p2 is their relative momentum. |0q(»")| 2 is 
a weight factor, which takes care of the indistinguisha- 
bility of identical particles. Sp(r) is the source function. 
The integral transform in Eq. (1) does not have a unique 
inverse due to the "on shell" constraint P ■ q = [3. 
Therefore, to extract the distribution Sp(r) from the ex- 
perimentally measurable correlation C(P,q) additional 
assumption must be made on the functional form of the 
source function Sp(r), which should be consistent with 
the symmetry of the system. 

Two kinds of symmetry are related to the present prob- 
lem. The first kind is the rotational symmetry, i.e. the 
symmetry property with respect to the space rotation. 
For a central heavy ion collision there is axial symmetry 
about the beam axis, defined as longitudinal direction. 
This amounts to isotropic in the transverse plane, i.e. 
the source appears as a circle in the transverse plane 
with the same radius in all directions. For a non-central 
or peripheral collision, the axial symmetry is violated. 
The source has an ellipse form in the transverse plane. 
The radii along the two main axes are nonequal. 



The second kind of symmetry is the inversion symme- 
try, i.e. the symmetry property with respect to the space 
inversion Xi — > —x%. This symmetry holds for all the 
three coordinates in the laboratory frame. However, in 
the present relativistic heavy ion experiments the HBT 
analysis is carried out not in the laboratory frame but 
in the Bertsch-Pratt frame 30 — Xi,x ,x s , where the 
longitudinal component xi still points along the beam di- 
rection; the outward component x Q is defined as parallel 
to the projection K of the pair momentum P on the 
transverse plane; the sideward component x s is perpen- 
dicular to both xi and x Q . This coordinate system has, in 
particular, the advantage that the ratio x /x s is sensitive 
to the duration of hadronization [j| 0] It is essentially 
a pair-by-pair rotation of the laboratory frame so that 
the new coordinate x points to the pair transverse mo- 
mentum K for all the pairs, cf. Fig. 1. 

In the Bertsch-Pratt frame, the inversion transform of 
the coordinate x„ 



(2) 



will change the direction of the vector K, and there- 
fore, is no longer a symmetry transform of the system. 
In other words, the existence of a fixed vector K along 
the direction x violates the inversion symmetry about 
x = 0, and a Gaussian source becomes invalid along 
this direction. 
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FIG. 1: The pair- by-pair rotation from labora- 
tory frame to Bertsch-Pratt frame. In the figure 
■di is the rotation angle of a;i a b for the «th pair. 
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FIG. 2: Upper row: the x -x s distribution of hadron 
pairs in a typical event of Au-Au central collision at 
^/snn = 200 GeV from AMPT simulation at three dif- 
ferent time. From left to right t = 5, 15, 25 fm, respec- 
tively. The lower row is the projection to the x Q axis. 



The aim of the present letter is to discuss how to take 
the inversion asymmetry along x into account in con- 
structing a source function. 

The inversion symmetry is an internal property of the 
distribution, which is defined in the intrinsic frame where 
the average value of coordinate vanishes, 



xP(x)dx = 0. 



(3) 



Thus by definition a Gaussian is always inversion sym- 
metric. 

The inversion-asymmetry degree A can be defined in 
the above-mentioned frame as 



A = 



J^lPW-Pj-xtfdx 
\ 2/_°° 00 [PW] 2 dx 



(4) 



which vanishes for inversion symmetric distribution and 
equals unity for maximum asymmetric distribution, 
where P(x) ^ only when P(—x) — 0. 

In current relativistic heavy ion experiments while per- 
forming HBT analysis a Gaussian source function is usu- 
ally used. Such a function may be isotropic or anisotropic 
consistent with the rotational symmetric or asymmetric 
of the collision, but could never be inversion asymmetric. 



As stated above, in relativistic heavy ion collisions the 
inversion symmetry holds for the source distribution in 
the laboratory frame, where the axes in the transverse 
plane are fixed and are common for all the particle pairs, 
but is broken along the x direction of the Bertsch-Pratt 
frame, cf. Fig. 1, so the Gaussian source function is no 
longer valid along this direction. 

In order to illustrate how the x distribution vio- 
lates the inversion symmetry and how the inversion- 
asymmetry develops in the heavy ion collision process, 
we take a Monte Carlo generator AMPT J\ as exam- 
ple. This model is based on non-equilibrium transport 
dynamics. It contains four main components: the ini- 
tial conditions, partonic interactions, conversion from 
the partonic to the hadronic matter and hadronic in- 
teractions. The initial conditions, which includes the 
spatial and momentum distributions of minijet partons 
from hard processes and strings from soft processes, are 
obtained from the HIJING model [8] in which eikonized 
parton model is employed. The time evolution of partons 
is then treated according to the ZPC |9| parton cascade 
model. After partons stop interacting, a combined coa- 
lescence and string fragmentation model is used for the 
hadronization of partons. Scattering among the resulting 
hadrons are described by a relativistic transport (ART) 
model [T(| which includes baryon-baryon, baryon-meson 
and meson-meson elastic and inelastic scattering. 

In the upper row of Fig 2 is sketched the x a -x s distri- 
bution of hadron pairs in a typical event of ^/snn = 200 
GeV Au-Au central collision at three different times - 
t — 5, 15, 25 fm. The inversion-asymmetry along x can 
clearly be seen. The asymmetry degree develops stronger 
and stronger with the increase of time. The lower row is 
the projection to the Xq axis. 

In order to extract the asymmetry information of 
source we need a source distribution function which can 
extrapolate from nearly symmetric, nearly Gaussian, to 
highly asymmetric. For this purpose we propose the fol- 
lowing three-parameter A, B, a distribution 



P(x) = < £ N+1 r(iv+i) ( a T x "> Ne 
|0 



±x < a, 
±x > a. 



(5) 



Transforming to the intrinsic frame, satisfying Eq. (3), 
we have 



where 



= fsgn(B) • '^g-g" for sgn(B) • x < sgn(B) • (N + l)B, 

m 1 for sgn(B) • x > sgn(B) ■ (N + l)B, 



N>0, sgn(B)= +j (7) 



r 



In Fig. 3 are plotted the P/vB a (aO's for three different N. It can be seen that the distribution varies from ap- 
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FIG. 5: (Color on line) The fit of the x a distri- 
FIG. 3: (Color on line) The NBa distributions for 3 dif- bution of the AMPT event at t = 25 fin used in 
ferent values of N. The absolute values of B are adjusted Fig. 2 to NBa (solid line) and Gaussian (dashed line), 
to make the 3 curves have successively lower height. 
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FIG. 4: The asymmetry degree of NBa distribution. 



TABLE I: The size R a and asymmetry degree 
A of a typical -^/snn = 200 GeV Au-Au cen- 
tral collision event at t = 5, 15, 25 fin used in 
Fig. 2, obtained from the fit of the x Q distribution 
to the NBa and Gaussian distributions, respectively 





Distrib 
NBa 


ution 
Gaussian 


Time 


Ro A 


Ro A 


5 fm 


2.28 0.0065 


2.28 - 


15 fm 


3.30 0.0565 


3.30 - 


25 fm 


4.90 0.1137 


4.90 - 



square of coordinate x is 



Rnbo, = 



v 2 PNBa(x)dx = \B\ ■ VN + 1. (9) 



proximately symmetric about x Q — at large N to highly 
asymmetric at small N. At sufficiently large N the NBa 
distribution Eq. (6) can mimic the Gaussian distribution 
very well. 

It turns out that the asymmetry degree A of this dis- 
tribution 



r(i)r(jv + i) 

Y{2N + l)T(N + §) 



[2(JV + l)]2JV+l e -2(iV+l) 
(8) 

depends only on the parameter N but not on B and a. 
The asymmetry degree A as function of N is plotted in 
Fig. 4. It can be seen that A varies from ~ at large N 
to ~ 0.5 at small N. 

The size, or "radius", R defined as the root mean 



We fit the x distributions of the Monte Carlo events 
shown in the lower row of Fig. 2 to the NBa distribu- 
tion (6) and present the results for t — 25 fm in Fig. 5 
as solid line. For comparison the same distribution is 
also fitted to Gaussian and shown in the same figure as 
dashed line. It is evident that the NBa fit can account 
for the inversion-asymmetry of x distribution present in 
the Monte Carlo data, while the Gaussian fit fails to do 
so. The size R and asymmetry degree A resulting from 
the fit are listed in Table I. Using the NBa source func- 
tion we are able to extract the asymmetry degree of the 
source, which could not be achieved using a Gaussian 
source function. As for the source size the results from 
both NBa and Gaussian source functions are the same. 

The inversion asymmetry exists only in the outward 
direction, so we assume Gaussian sources in xi and x s , 
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while NBa source in x Q , 



Performing Fourier transform we get 



3 -|(g?-R?+«Jfl;) e -ig„(iV+i)B 



(1 - iq B) N +^ 
The correlation function is 
C 2 (q) = 1 + A|p(q)| 2 
= 1 + A 



(10) 



(11) 



[1 + q*B x \V f + 1 ) 



-(9i 



(12) 



Through fitting the experimentally measured correla- 
tion function to Eq. (12) the parameters Ri, R s , N and i? 
are obtained. From these parameters beside being able 
to extract the size parameters Ri and R s in the longi- 
tudinal and sideward directions, we are also armed with 
a powerful tool for extracting the asymmetry degree A 
and size parameter R a in the outward direction by using 
Eq's. (8) and (9). 

It is argued that under the Bertsch-Pratt conven- 
tion 00 popularly used in the HBT analysis of rela- 
tivistic heavy ion collision there is a fixed vector — the 
total transverse momentum K of the pair, along the out- 
ward direction. The existence of such a vector violates 
the inversion symmetry along this direction, making the 
conventional Gaussian source function invalid. Although 
the size of source can be extracted from the experimen- 
tal data using a Gaussian source function, an important 



feature of the system 
is missed. 



the asymmetry degree along x a 



In the present letter the inversion-asymmetry along 
x Q is pointed out in the first time. A 3-parameter 
source distribution - - the NBa distribution is pro- 
posed, which can extrapolate from nearly inversion- 
symmetric to highly inversion-asymmetric. Using a 
source function, which is Gaussian in the longitudinal 
and sideward directions, while NBa in the outward di- 
rection, the asymmetry degree A along the outward di- 
rection can be extracted from the experimental data to- 
gether with the source-size parameters i?;, R s and R D . 

The inversion asymmetry along x in a typical Au- 
Au central collision event under the Bertsch-Pratt con- 
vention has been checked using the AMPT Monte Carlo 
model, and the successful extraction of the asymmetry 
degree A from the data is exhibited. To apply the pro- 
posed method to real experimental data is highly recom- 
mended. 
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